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ABSTRACT 

We develop two methods to estimate the bulk Lorentz factor of X-ray flare outflow. In the first 
method the outflow is assumed to be baryonic and is accelerated by the thermal pressure, for which 
the final bulk Lorentz factor is limited by the outflow luminosity as well as the initial radius of the 
outflow getting accelerated. Such a method may be applied to a considerable fraction of flares. The 
second method, based on the curvature effect interpretation of the quick decline of the flare, can give 
a tightly constrained estimate of the bulk Lorentz factor but can only be applied to a few giant flares. 
The results obtained in these two different ways are consistent with each other. The obtained bulk 
Lorentz factor (or just upper limit) of the X-ray flare outflows, ranging from ten to a few hundred, is 
generally smaller than that of the Gamma-ray Burst (GRB) outflows. 

Subject headings: Gamma Rays: bursts— radiation mechanisms: nonthermal— X-rays: general 
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1. INTRODUCTION 

Cosmic Gamma-ray Bursts (GRBs) are the most lu- 
minous events ever known in the universe after the Big 
Bang. The Luminosity, duration a nd bulk Lore ntz fac- 
tor are crucial parameters of GRBs ([Piran 1119990 . Based 
on the duration GRBs can be divided into two main 
groups, short and long bursts. Once the redshift of a 
burst is measured, with a given cosmology model the 
luminosity distance is known, so are its luminosity and 
isotropic-equivalent 7— ray energy. Though the physi- 
cal processes producing prompt emission and afterglow 
of GRBs are still unclear, it is widely accepted that the 
outflows move ultra-relativistically (i.e., its velocity is 
very close to the speed of light). The most reliable evi- 
dence is the measurement of superluminal movement of 
the r adio afterglow image of GRB 030329 ()Tavlor et al. I 
2004). In the literature, several methods have been de- 
veloped to constrain the initial bulk Lorentz factors of 
GRB outflows. (1) The detection of the high energy 
photons imply the emitting region is optical thin for pair 
annihilation. T his argument leads to a lower limit on 
Loren tz factor (|Shemi fc Piran I IT990t ILithwick fe Saril 
2001). (2) In some cases the peak in the early GRB af- 
terglow lightcurves marks the deceleration of the external 
forward shock, at which the Lorentz factor is about half 
of the initial value. In such cases the initial Lorentz fac- 
tor, weakly depending on the isotropic-equivalent kinetic 
energy a nd the density of the medium, can be tightly con- 
strained |SirrSEran][i9p ^^^T^^T^M) ■ (3) 
Detection of a thermal emission component in the GRB 
afterglow spectrum provides a relatively direct way to es- 
timate the Lorentz factor. Assuming a thermal radiation 
efficiency, with the measured temperature and flux of the 
therma l component , one can calculate the bulk Lorentz 
factor (|Pe'er et al. 1 120071) . (4) The non-detection of a 
hard X-ray to soft gamma-ray background emission in 
GRB prompt stage gives an up per limit on the in itial 
Lorentz factor of GRB outflows |Zou fe PirarHIMpf) . 

Many Swift GRBs were followed by energetic X-ray 



yzfanOpmo.ac.cn(YZF), dmwei@pmo.ac.cn(DMW) 

1 Purple Mountain Observatory, Chinese Academy of Sciences 



flares. Usually their fluences are only about 1 — 10% 
of that of prompt emission. However in quite a few 
events, the energy of X-ray flare is comparable to that 
of prompt emission. The temporal behavior and the 
hardness ratio evolution of X- ray flares are similar t o 
those of prompt emission pulses (jChincarini et al. I[2010l) . 
supporting the idea that they have the same physical 
origin of the prompt emission, i.e., they are also due 
to the activity of c entral engine (e.g. iFan fe Weil [20051: 
IZhang et al. | |2006). In such a kind of model, the cen- 
tral engine launches new outflow at late times. As in 
the prompt emission phase, the bulk Lorentz factor of 
the flare outflow is a crucial parameter. Unfortunately, 
most models constraining the bulk Lorentz factor of the 
GRB outflow are invalid for t he X-ray flares. For e xam- 
ple, the methods de veloped in Sa ri fc Piran ] (fl999h and 
IZou fc Piranl (I20TI are irrelevant since there is a pre- 
ceding and more energetic GRB outflow expanding into 
the medium. So far there are three kinds of specula- 
tions on the bulk Lorentz factor of the flare outflows: (a) 
The typical bulk Lorentz factor of flares is just tens and 
is considerably small er than that of th e outflow power- 
ing prompt emission (|Fan fc Wei 1120050 : (b) The typical 
bulk Lorentz factor of flares is higher but not much higher 
than that of the outflow giving rise to pr ompt emission 
(|Burrows et al. 1 120051: IZhang et all 120060 ■ Please note 
that in both case (a) and case (b), the flare photons 
are powered by the energy dissipation within the newly 
launched outflow, (c) In the X-ray flare model of up- 
scattered forward shock emission, l ate outflow with a 
bulk Lorentz factor ~ 10 5 is required (Panaitcscu 2008). 
The divergency between these arguments are very large. 
In this work, we develop two different methods, as de- 
scribed in section [21 to estimate the bulk Lorentz factor 
of the outflows. The case studies are presented in sec- 
tion [3] We summarize our result with some discussion in 
section 2] 

2. THE METHODS 

Method I: The physical composi tion of X-ray flare 
outflows is not well constrained yet. IFan et aT~l (2005) 
suggested that in the neutron star— neutron star (or black 
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hole) merger model the outflow powering X-ray flares fol- 
lowing short GRBs might be Poynting-flux dominated 
since the fallback accretion onto the newly formed black 
hole is too low to launch an energetic ejecta via neu- 
trino annihilatio n. But for long GRBs t he argument is 
weak. Recently iChincarini et al.l ([2010L and the refer- 
ences therein) argued that the X-ray flare outflow was 
likely Poynting-flux dominated if the prompt emission 
was powered by the magnetic energy dissipation. How- 
ever the physical origin of the prompt emission is still in 
debate. Therefore here we assume that the flare outflow 
is baryonic. The central engine releases a large amount 
of energy into a compact region. Therefore the outflow 
is very hot with a temperature ~ 1 MeV, depending on 
the total luminosity of the flare outflow (L) as well as 
the initial radius of the outflow getting launched (Rq). 
As in the case of fireball powering prompt emission, the 
flare outflow will be accelerated by the thermal pressure 
until it becomes optically thin (i.e., the optical depth 
r ~ 1) or saturates at a radius i?f, depending on the 
outflow is baryon-rich or not. In the baryon-poor case 
the outflow becomes transparent at Rf < t/Rq, where rj 
is the dimensionless entropy of the initial ejecta. The 
thermal energy has not been effectively converted into 
the kinetic energy of the outflow and will give rise to a 
quasi-thermal radiation component. The absence of such 
a soft component in the data may suggest that the flare 
outflow is baryon-rich, for which the final bulk Lorentz 
facto r can be estimated as T x ~ 77 ~ Rf/Ra (Piran 
1999). On the other hand, the optical depth of the pho- 
ton at the radius Rf can be estim ated as (|Paczvhski I 
1199(3: IDaiene fc Mochkovltch1[200l : 

POO 

t ~ / (1 - (3 x )n<r T dR ~ 1, (1) 

JRt 

where n ~ L/4irR 2 T x m p c 3 is the number density of elec- 
trons coupled with protons in the observer's frame, <tt 
is the Thompson cross section, and m p is the rest mass 
of protons. Combing with the relations Rf ~ T x i?o and 
/3 X ~ 1 - l/2r x , eq.® gives 



8TrT^.m p c 3 Rq 

Therefore the bulk Lorentz f actor of X-ray flares is 
related to L and Rq as (e.g . iMeszaros fc Reesl 120001 : 
INakar et al. 1120051: lFa^~|[20Ta) 

Tx < Tma X = 5 X 10 2 L 5 g Rq q , (3) 

throughout this work, the convention Q x = Q/10 x has 
been taken into account except for some special nota- 
tions. For rj <C r max , the photospheric radius of the 
flare outflow is estimated to be ~ 6x 10 10 cm L^-q^ 3 
(jPaczvhski I I1990L iFan 1 120101 ) . much larger than R { ~ 
10 s cm 772-^0,6- With a proper R and the observed X- 
ray flare luminosity L x (<~ e x L, where e x is the X-ray 
flare efficiency and is taken to be compar able to that of 
GRBs, i.e., - 0.1 (e.g. lFan fc Rrar7l[200l ). we are then 
able to give an upper limit on T x . The advantage of this 
method is that it may apply to a good fraction of X-ray 
flares. The limit of this approach is that it is only valid 
for the baryonic outflows. 



Method II: Th e quick decline o f the X-ray flares 
(|Piro et al. I 120051 iBurrows et al~l [2001 may have 
impo sed a tight c onstr ai nt on the emission radius 
R^ flZhang et al. I 120061 lLazzati fe Begelman I 120061 : 
iDai et al. I 12007ft provided that the quickly decaying 
X-ray emis sion is the high latitude component of the 
flare pulses ([Fenimore et al. 1996; K umar fe Panaitescu I 
2000, in some literature it is called the "curvature ef- 
fect"). This interpretation, i f correct, require s a very 
large variability timescale ST ([Fan fe Wei 1120051 see also 
below). The corresponding bulk Lorentz factor thus can 
be estimated by T x w [R x /{2c5T)] 1 / 2 , please note that 
the timescales involved here and below are all in the rest 
frame of the GRB source. If R x as well as ST can be 
reliably estimated, T x can be reasonably constrained. 

For simplicity, we take the leading late internal shock 
model for illustration. Assuming that there are two 
shells, the slow one is ejected at T and moves with a bulk 
Lorentz factor r s , the faste r one is with Tf, and thei r ejec- 
tion interval is St (see also fkobavas hi et al. I IT997). The 
widths of both shells are ~ A (in the observer's frame). 
The faster one would catch up with the slower at a radius 

flcon « 2T 2 s cSt. (4) 

After the merger, the newly formed shell has a bulk 
Lorentz factor T x ~ yTf r s ~ if the mass of these two 
shells are comparable ([Piran 1 11999ft . The timescale of 
that merger is determined by the reverse shock crossing 
the fast shell and is ~ A/c. The corresponding radius is 

R x ~ (2r 2 s c5t + 2rlA), (5) 

at which the emission peaks. In the internal shoc k phase, 
usually the electrons cool rapidly ([Piran 1 11999ft . After 
the ceasing of the internal shock, the em ission is from the 
high l atitude 9 > 1 /T x and takes a form dFenimore et al. I 
[19961 iKuniar fe Panaitescu II2IM IFan fe Wei 112005^ 

Foe (T/6T)-( 2+ P\ (6) 

where ST « R x /(2Tlc) « A/c for T 2 > T 2 and 
St ~ A/c, and the standard convention of the flux as 
a function of both time and frequency f v oc T~ a v~fi 
is adopted. In this work, we assume that the ejecta is 
uniform. 

In reality, the flare consists of many (for example k) 
pulses. After the ceasing of internal shocks at T t , what we 
observe is the high latitude emission of the early pulses. 
The X-ray flux declines as 

k 

F X w X>*,<[(T - ToJ/STi]- {2+f> -'\ (7) 
i=l 

where i represents the i— th pulse, F Vxt % is the peak emis- 
sion of the i— th pulse, Tqj is the time when the i— th 
pulse is ejected. Such a decline is much steeper than 
(j7j;)-(2+/3 ) as i ong as T% > ma x{STi} (see Fig.2 of 
IFan fc Wei I ([2005) for illustration). This can be under- 
stood as follows. In the popular internal shock scenario, 

2 Please bear in mind that the following discussion is valid as long 
as the prompt emission consists of many pulses. These pulses could 
be powered by either late internal shocks or late internal magnetic 
energy dissipation. So the validity of our method II, different from 
method I, is independent of the physical composition of the outflow. 
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each prompt pulse is independent and is emitted at a 
radius ~ i? x . At late times, the emission contributed by 
early shells is from a very large angle, and is very weak 
due to the relativistic beam effect. The observed flux 
is thus dominated by the curvature emission component 
of the last pulse (if all STi are comparable) or one early 
pulse having a very large duration ~ T t . As a result, the 
net flux of these shells can be approximated by 



Fx 



F^, k [(T-T 0ik )/ST k ]-^\ 



"(2+/3) 



(8) 



since 



which could be far more sharply than (T/T t ) 
the internal sho cks cease- 
As shown in iLiang et al. I (|2006[ ) and iZhang et al. I 
(|2007l) . the sharp decline of some X-ray flares are 
well fitted by T x cx (T - T n )~ ( 2+/3) , which is consis- 
tent with the curvature ef fect (|Fenimore et al. 1 119961 : 
iKumar &: Panaitescu I [2000) interpretation of that phe- 
nomenon. The inferred To is a good approximation of 
the ejection time of the last dominant pulse (rather than 
the re-activity time of the central engine unless there is 
just a pulse). 

As shown in eq.®, for the last dominant pulse, ST ~ 
T p — T), where T p is the peak time of the last dominated 
pulse. On the other hand, the curvature emission com- 
ponent (i.e., the tail emission) lasts T ta ii until the edge 
of the ejecta is visible (see FigJ3]). For an ejecta hav- 
ing a half-opening angle f5j (which can be estimated from 
the achromatic breaks of the late afterglow lightcurves 
(jRhoads l fl999^l . the prompt emission radius can be es- 
timated as 

i? x ~ _ g^L_ (9) 



So we have 



Ttail 



(l-costfj)- 



11/2 ^ r ftail 11/2 



L 2(l-COS#j)(ST J 



L 2(l-cos0j)(Tp-To) J 



( 10 ) 

In many cases, it is inconvenient to get T p and To, respec- 
tively. Fortunately, a simpler way to estimate T ta ii/<5T 
is available. We introduce the factor 1Z to denote the 
ratio between the peak flux of the X-ray flare (F XtP ) and 
the flux when a cutoff emerges (T XjC ), i.e., 1Z = T XjP /F XiC 
(see Fig El). With eq.©, we have ' 

( Tail 2+ ^ Tail^ 

Combing eq. (|TU]) with eq. ([TTl) . we have 



(12) 



With a typical £ ~ 1, T x oc ft 1 / 6 . So our estimate of T x 
can not be modified significantly by the uncertainty of 
TZ. The choice of 9y instead, is crucial for our purpose. 
For the X-ray flare outflow, there is no simple/reliable 
method to estimate its half-opening angle f3j fl aro . How- 
ever, it is reasonable to argue that (9j fl arc can not be 
much smaller than 6\grb. This is because X-ray flare(s) 
have been observed in about 40% Swift GRBs, so sta- 
tistically speaking, 0j jflare ~ n/oTs^grb ~ 0.76> j:GRB (X. 
Y. Dai 2006, private communication), where 0j,GRB is 
the half-opening angle of the GRB ejecta that can be 
inferred from the late afterglow break. Such a small 
correction can just increase our estimate of T x (see 
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Fig. 1. — Upper limits of the bulk Lorentz factor of outflow pow- 
ering X-ray flare in GRB afterglow and the center time of the 
X-ray flares. The two types of triangles represent the upper limits 
on the bulk Lorentz factor estimated in method I. The hallow and 
filled triangles are for the single power-law spectrum assumption 
and the Band function spectrum assumption, respectively. Stars 
arc the bulk Lorentz factor, estimated in method II, of two giant 
flares. 



eq.(|T2j) for details) by a factor of 1.4 and thus does 
not change our basic conclusion. As a consequence, 
in the case studies (see below), we simply take 6*j = 
^j.flare = #j,GRB- The error of the resulting T x can be 
estimated straightforwardly and can be approximated by 

ST X ~ r xV / 3(5 ^ F (^) 2 + { 5 4f ) 2 , where STZ and 66, are 

the errors of TZ and 6j, respectively. For (3 ~ 1 and 
51Z/1Z < 6<5fV0j, we have 5Y* ~ Y x 59j/9y For a conser- 
vative estimate 89J9- } ~ 1/2, we have <5r x ~ 0.5r x . 

Different from method I, the current method can give 
a somewhat reliable estimate of T x rather than an upper 
limit. However, a cutoff at the end of the quick decline 
is needed to achieve that goal. Unfortunately an unam- 
biguous cutoff in the quick decline phase has just been 
identified in a few giant flares, limiting the application 
of method II. We'd like to point out that method II is 
also valid to estimate the bulk Lorentz factor of the late 
outflow of GRBs if the quick decline phase of prompt 
emission is also attributed to the so-called curvature ef- 
fect. 

3. CASE STUDIES 

For method I, we take a sample consisting of 36 X-ray 
flares detec ted in 14 GRBs with m easur ed redshift, as 
report ed in iChincarini et al~l (j2007j ) and iFalcone et al.~l 
(2007). The average bolometric luminosity, i.e., the bolo- 
metric energy divided by the duration, is derived for both 
the power-law and t he Band function fits to the flares 
(|Falcone et al. II2007D . As already mentioned, the total 
luminosity of the flare is taken to be 10 times than that of 
the X-ray emission. Rq is taken as 10 7 cm, which is com- 
parable to the radius of a neutron star or the last stable 
circle orbit for a rapidly rotating black hole (In the ther- 
mal radiation modeling of some GR Bs, Rp is estimate d 
to be ~ 10 cm or even larger (e.g. iPe'er et al. 1 120071) . 
So our estimate is likely conservative). The results are 
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Fig. 2. — The distribution of the upper limits of bulk Lorentz 
factor of outflow powering X-ray flare in GRB afterglow. Two types 
of slant lines are for the single power-law spectrum assumption and 
the Band function spectrum assumption, respectively. 



presented in Fig[T] and Figj^J The resulting upper limits 
on T x are between tens to hundreds, which are gener- 
ally smaller than the initi al bulk Lorentz factor of GRB 
outflows ~ 10 2 - 10 3 r e.g.lPiranlll99<llLithwick fc"SarT1 
I200H IMolinari et aT~l 120071 IZou fc Piran 1 12010^ For 
the baryonic outflow, we have T x ~ L/Mc 2 , where 
M ~ 5 X 10 _6 i5o T~\Mq s _1 is the mass loading rate of 
the outflow. The inferred T x > 10 following the flares in 
GRB 050502B and GRB 050724 suggests that the flare 
outflows are still relativistic, strengthening the connec- 
tion between the X-ray flares and the prompt soft 7— ray 
emission. Since the typical total luminosity of the flare 
outflows is about two or more orders of magnitude lower 
than that of the GRB outflows, the mass loading rate 
of the flare outflow is expected to be lower than that of 
the GRB outflow otherwise relativistic ejecta can not be 
launched. This is reasonable since the pole region of a 
dying massive star (or the remnant of the merger of two 
compact objects) is expected to be more and more clean 
as the time goes on. 

Below we focus on method II. The sharp decline has 
been well detected in a good fraction of X-ray flares. 
For our purpose, a cutoff, emerging when the emission 
at the edge of the ejecta is within the line of sight, 
of the tail emission is needed to constrain -Rfiarc reli- 
ably. Such a cutoff has just been possibly identified in a 
few cases, for examp le, the giant flare in GRB 050502B 
([Falcone et al. 112006ft. a nd the giant flare in GRB 050724 
( Campana et al. 1 12006). The reason for such a rare de- 
tection is the following. In reality, the X-ray emission is 
not only contributed by the high latitude emission but 
also contributed by the forward shock emission. The 
latter would dominate over the former when the high 
latitude emission component has dropped by 2 or more 
orders (correspondingly, T ta y ~ 5ST for j3 ~ 1) unless 
the X-ray flare is strong enough, or the forward shock 
emission is very dim. On the other hand, when the flux 
has dropped to a level of ~ 10~ n erg s _1 cm~ 2 , the sig- 



FlG. 3. — A schematic plot of the curvature emission component 
(tail emission) that is useable for method II. 

nal/noise ratio is not high enough to get a good quality 
detection and thus renders the identification of the cutoff 
difficultly. The data (together with corresponding refer- 
ence) are presented in the Table fT] 

Note that there is a possible independent argument fa- 
voring aT x ~ tens for giant flares f ollowing GRB 0505 02b 
and GRB 050724. As suggested in lFan fc Weil (|2005ft . in 
the leading late internal shock model, the X-ray flare out- 
flow with a r x ~ tens will catch up with the initial GRB 
outflow when the latter has swept a large amount of ma- 
terial and then got decelerated. Such an energy injection 
process would give rise to a flattening (if there is also a 
wide range of the bulk Lorentz factors of the X-ray flare 
outflow) or re-brightening signature (if the range of the 
bulk Lorentz factors of the flare outflow is narrow). The 
re-brightening has been well detected in the late X-ray af- 
terglow of GRB 050502b and GRB 050724. To attribute 
the very late X-ray flare in GRB 050502b to an energy in- 
jection ca used by the giant X-r ay flare outflow, a T x ~ 20 
is needed (|Falcone et al. If2006h . matching our re sult per- 
fectly. As for GRB 050724. IPanaitescu I (poM ) showed 
that the late multi-wavelength re-brightening could be 
well reproduced by an energy injection. Before the en- 
ergy injection, the kinetic energy of the GRB ejecta 
Ek ~ 10 50 erg and the number density of the medium 
is n ~ 0.1 cm~ 3 . The decelerating GRB ejecta has a 
bulk Lorentz factor - 10 (£ k:5 o/n-i) 1/8 (7yi0 4 s)- 3 / 8 . If 
this energy injection is also caused by the flare outflow 
catching up with the GRB ejecta, aT x ~ 10 is needed, 
which is consistent with our result. 

Finally, we would like to point out that only the X-ray 
flares brighter than the forward shock X-ray emission 
can be identified. There could be some faint X-ray flares 
which have even lower bulk Lorentz factor. 

4. DISCUSSION AND SUMMARY 

Bright X-ray flares have been well detected in the af- 
terglow of a considerable fraction of GRBs. The radia- 
tion mechanism powering these events is still unclear. A 
widely accepted hypothesis is that new outflow should be 
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TABLE 1 

The physical parameters, including the estimated bulk Lorentz factor, of giant flares detected in GRB 050502b and GRB 

050724. 



GRB tc 
050502b 152s 
050724 122s 



n 



r x (i) r x (ii) 



Reference 



1.3 
1.34 



1100 
72 



■ 0.13 
< 0.15 



< 666 

< 120 



22 
13 



launched by the central engine and the flares should be 
powered by the energy dissipation within the new out- 
flow. The lack of an optical flare associated with the 
X-ray flare is in g eneral consistent with such a scenario 
([Fan fc Wei 1120051) . The bulk Lorentz factor of the newly 
launched outflow is a crucial parameter for us to under- 
stand X-ray flares. For example, with the inferred bulk 
lorentz factor we can see whether the outflow is relativis- 
tic or not. Then we can estimate the corresponding mass 
loading rate in the pole region of the central engine. If the 
inferred bulk Lorentz factor is so high that surpasses the 
upper limit given in eq.([3]), the baryonic outflow model 
will be ruled out. 

In this work we have developed two methods to esti- 
mate the bulk Lorentz factor of X-ray flare outflow (see 
section 2). In the first method the outflow is assumed 
to be baryonic and is accelerated by the thermal pres- 
sure, for which the final bulk Lorentz factor is limited 
by the outflow luminosity as well as the initial radius 
of the outflow getting accelerated. Such a method may 
be applied to a considerable fraction of flares. How- 
ever it can only give an upper limit and is invalid if 
the flare outflow is Poynting-flux dominated. The sec- 
ond method, based on the curvature effect interpretation 
of the quick decline of flare, is independent of the phys- 
ical composition of the outflow and can give a better 
constrained estimate of the bulk Lorentz factor of the 
flare outflow but can only be applied to a few events 0- 
The obtained bulk Lorentz factors (or just upper limit) 
of the X-ray flare outflows in these two different ways 
are consistent with each other and are generally smaller 
than that of the GRB outflows (see section 3). How- 
ever, the flare outflows are still relativistic and the corre- 
sponding baryon loading rate is expected to be very low 
(M < 10~ 5 Mq s" 1 ). This finding, together with other 
results from the analysis of the X-ray flare data such as 



Falcone et al. (2006) 
Campana ct al. (2006); Panaitcscu (2006) 



the temporal behavior and the hardness ratio evolution, 
strongly favors the hypothesis that X-ray flares are the 
low energy analogy of the prompt so f t 7— ray emission 
dFan fc Wei |[200a IZhang et al. Il200l IChincarini et al. I 
12010b IMareutti et al. I I2010t IShao et al. 1 120101) . With 
method II, the Lorentz factors of the two giant flares 
are found to be ~ a few x 10, well below the upper limits 
given by method I (see also Tab.l). Hence the baryonic 
outflow model is not challenged. For the magnetic out- 
flow with reasonable baryon loading, a moderate Lorentz 
factor is also possible. Further data like the high linear- 
polarization of the X-ray photons is needed to claim the 
magnetic nature of the flare outflow. 

If the X-ray flare outflows do have a bulk Lorentz factor 
~ 10 5 , as suggested i n the flare model o f up-scattered for- 
ward shock emission (jPanaitescu 1120081) , the acceleration 
can not be due to the thermal pressure. The only known 
such a kind of extreme astrophysical object is the pulsar 
wind, which may move with a Lorentz factor as high as 
~ 10 6 and its acceleration is likely a result of the signif- 
icant magnetic energy reconnection (|Kennel fc Coroniti I 
Il984f ). It's less likely to be the case for the X-ray flares 
since the central engine of GRBs is not in a cavity as 
clean as that of a pulsar at an age of > 10 3 years. As 
for the flare model of up-scattered forward shock emis- 
sion, one more challenge is that the similarities between 
the GRB prompt emission and the X-ray flare emission 
strongly suggest a similar physical origin of these two 
kinds of phenomena. 
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